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. We apply an improved renormalization group analysis for pure Yang-Mills theory in 

' four dimensions in the Landau gauge and 0{N) nonlinear sigma model in two dimensions. 

' We find a pole in the gluon two-point function at one-loop order, which is generated non- 

, perturbatively. In Yang-Mills theory this may correspond to the expectation value of a gauge 

I . invariant operator and emerges as Mp/Aj^ ~ 0.81, where A^fg is the MS scale. 

, §1. Introduction 

■ The improved perturbation theory formulated by Geh-Mann and Low^^ is a the- 
, oretical framework with which, using the ideas of the renormahzation group with the 
' results of perturbation theory to a given order, one can determine something about 

■ the next order of perturbation theory. Employing the same type of philosophy, vari- 
i ational approximation methods, by which one can extract certain nonperturbative 
' information from the results of perturbation theory, have been developed in both 

theoretical and numerical approaches. 

D i The essential idea of the variational approximation scheme, which is often called 

' an improved mean field approximation in the literature, can be summarised simply 

■ by the expression "the principle of minimal sensitivity" .^^ Suppose that we have a 

' I A 

^ , Lagrangian £ of interest. For concreteness, we use the (j) theory as an example. 

\—{ ' We first introduce a mean-field Lagrangian Cm = m'^cjp' /2, where is a parameter 

to be tuned such that the quantity we would like to evaluate is independent of this 
parameter. We then rewrite the original Lagrangian as 



where we have introduced the formal coupling constant A. We then regard the first 
two terms as the unperturbed action and the terms in the parentheses as a pertur- 
bation with respect to A, and calculate the quantity of interest using perturbation 
theory. This A will be set to 1 at the end of calculation, and thus, in principle, any 
physical quantity calculated with this Lagrangian cannot depend on m^, but, as we 
will see, the dependence remains if terms of higher order in A are dropped. Fi- 
nally, we tune the parameter for the calculated quantity to be independent of it 



typeset using VTpT^.ds (Ver.0.9) 
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as should be the case. In other words, one should choose the parameter for the quan- 
tity of interest to be on a plateau whose emergence signals that the approximation 
scheme works well. This is the meaning of the expression "the minimal sensitivity" . 

When we apply the variational approximation for a massless theory, it is con- 
venient to start with the massive counterpart £massive('TT'^i 9) to the massless theory 
and employ the ordinary perturbation theory with respect to the coupling g. All we 
have to do is to replace the mass and coupling as ^ (1 — A)m^ and g ^ Xg in 
the result of the perturbation theory and then keep the desired order in the formal 
coupling constant A and set A = 1. Here, the mass plays the role of the parameter 
to be tuned. 

There have been many works carried out with this method. For the case of 
anharmonic oscillators, it was shown that the new perturbation series obtained with 
this recipe is convergent. (See also Ref. 4).) Some progress has been made in 
the study of matrix models and reduced Yang- Mills models. ^^"^^-^ In Ref. 16), D- 
dimensional pure Yang- Mills theory is studied in the context of the reduced model 
and convergence of Monte Carlo data is demonstrated. In Ref. 7) and, subsequently, 
Ref. 8), the IIB matrix model, which is defined as the reduced maximally supersym- 
metric Yang-Mills theory and is also conjectured to be a nonperturbative definition 
of superstring theory, is studied, and the emergence of four-dimensional space-time 
as a vacuum of string theory is sugg estedEll 

Regarding applications to quantum field theory, in Ref. 17) the variational 
method is combined with the renormalization group method. The mass gap of the 
Gross-Neveu model in two dimensions is calculated and is found to exhibit good 
agreement with the exact solution. Furthermore, in a subsequent series of papers, 
the method is applied to QCD in order to study, among other things, the dynam- 
ical origin of the quark mass. In this paper, we call the method developed there 
"the improved renormalization group analysis" , and we review it in § 12.11 with some 
refinement. 

This improved mean field approximation can be regarded as an "improved Tay- 
lor expansion", which is a general scheme to improve the convergence of a Taylor 
expansion series of interest.®^ We note here that when we attempt to apply the 
improved mean field approximation to quantum field theories, we usually face the 
problem of determining at which stage we should make the replacement, i.e., how 
to make it compatible with the renormalization procedure. This problem arises be- 
cause both involve the choice of the unperturbed part of the action, and usually 
they do not seem to be compatible. The proposal made in 17) and 18) is that the 
replacement is carried out for the bare Lagrangian, as the original mean field ap- 
proximation suggests, and then the same renormalization scheme as for the original 
Lagrangian is applied. It was shown that this leads a consistent improved series; 
that is, the divergences are removed by this renormalization. Here we employ an- 
other approach. Inspired by the spirit of the improved Taylor expansion, we first 
apply usual renormalization scheme and improve the quantity of interest by using a 



*' In Ref. 12), higher-order calculations are given and the result strengthens the original idea. 
(See also Ref. 14) and references therein.) 
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renormalization group technique, and then apply the improved Taylor expansion to 
functions that admit a Taylor expansion. Surprisingly, our method gives the same 
improved function as the method in Refs. 17) and 18). We therefore believe that our 
general scheme still can be interpreted as an improved mean field approximation. 

In this paper, we apply the improved renormalization group analysis to the pure 
Yang-Mills theory in four dimensions in order to explore the nonperturbative na- 
ture of the theory, such as confinement. Though lattice theory provides a method 
for nonperturbative study, there is still a great need for analytic methods that we 
could trust. Recently, there has been growing interest in the condensate of a mass 
dimension two operator in QCD and Yang-Mills theory. Such condensation is be- 
lieved to be a key to understanding the confinement problem. I*)] The operator 
considered here is Z\ = ^ (volume) "^(minjy J d'^x{A^y), where U represents a gauge 
transformation, (volume) denotes the space-time volume, and "min" represents the 
operation of taking minimum with respect to gauge transformations. This operator 
is indeed gauge invariant, but non-local in a general gauge. In the Landau gauge, 
this is the expectation value of the gauge-variant local operator J (^^Afj_A^y^'^^ Then, 
in the Landau gauge, the mass term of the gluons, Tr m ^AnA'^, might be induced 



nonperturbatively through the four-point interaction, til This "induced mass" may 
appear as a pole of the gluon propagator. Thus, we perform the improved renor- 
malization group analysis to evaluate the mass pole of the gluon propagator in the 
Landau gauge. We do not demand that this "mass" be the physical mass of the 
gluon, since the mass term is not gauge invariant. However, what we calculate is as- 
sumed to be the expectation value of a gauge invariant quantity. Therefore, though 
this does not have the meaning of a mass, except in the Landau gauge, the quantity 
itself is thought to have a gauge invariant meaning. In other words, we simply pro- 
pose a possible way to calculate A in Yang-Mills theory. The condensation of these 
kinds of operators has been discussed by means of lattice simulations^^^ as well as 
from a phenomenological point of view.^^^'^-^^ 

This paper is organized as follows. In the next section we explain the improved 
renormalization group method developed in Refs. 17) and 18), with our own re- 
finement, employing the idea of the improved Taylor expansion. Then, in § [3] we 
demonstrate an application to the 0{N) nonlinear sigma model in two dimensions. 
In § m we calculate the mass pole of the gluon propagator in the Landau gauge by 
using the Curci-Ferrari model^^^ as a massive counterpart of Yang-Mills theory. We 
perform the calculation in a general covariant gauge, in which we can obtain the 
result in the Landau gauge by setting the gauge parameter ^ to zero. Section [5] is 
devoted to conclusions and discussion. 



*' For recent developments, see, e.g., Refs. 23), 24) and references therein. 

**' Note that no clear explanation has been given of the connection between mass generation in 
gauge theory and the condensate of the operators. There have been several arguments regarding 
this issue. (See Refs. 20) and 25).) Here we simply assume the mass generation of the Landau 
gauge. 
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§2. Improved renormalization group analysis 

We first describe an improved renormalization group analysis that we subse- 
quently apply to Yang- Mills theory on the basis of Refs. 17) and 18). We start with 
a Lagrangian density £(m^,g,/x), where is the mass parameter, 5 is a dimen- 
sionless coupling constant, and is a mass scale which is introduced to keep the 
coupling g dimensionless. This massive theory should be regarded as the massive 
counterpart of a massless theory of interest. As we have seen in the introduction, 
in order to obtain an improved perturbation series in A for a given quantity, such 
as propagators, in the massless theory, it is sufficient to replace m? with (1 — X)m? 
and g with A^f in the quantity calculated using perturbation theory with the massive 
action. After setting A = 1, one can, at least naively, expect that this new quantity 
describes that of the massless theory. Therefore, we first review the method for ob- 
taining a pole mass, in which we are interested here, improved by a renormalization 
group equation in the massive theory. As stated in the introduction, we realize this 
substitution by use of the improved Taylor expansion, which is in general expected 
to improve the convergence of a series and to extract information regarding the 
original function from its perturbation series of finite degree. We can easily apply 
this method to a perturbation series in quantum field theory, even after we have 
improved the perturbation series by use of the renormalization group technique. In 
this case, we simply improve the quantity of interest, which is written as a series in 
a coupling constant, by replacing the original parameters with new ones, including 
a formal expansion parameter denoted by A in the introduction. Here we note that 
there can be some quantities that do not admit a Taylor expansion, for example, 
a function of a fractional power in the coupling constant appearing through a non- 
perturbative refinement like the renormalization group. We do not consider such 
quantities and will improve only the parts that can be written as a power series 
with integer powers in the coupling constant. This is the general criterion for the 
improvement we propose in this paper. This procedure provides the same formula 
as that proposed in Refs. 17) and 18), where the improve mean-field approximation 
is applied to the bare Lagrangian, and then the renormalization and the renormal- 
ization group method are applied carefully. We thus conclude that our method is 
actually an improved approximation of a massless theory of interest. 

Here, let us start by explaining how to improve the pole mass expression in 
which we are interested in this paper. As explained in Appendix[Bl a pole mass can 
be written in terms of these parameters as 

n—L 

g'^ (2-1) 

where the quantities A^^l are constant, that is, are independent of m^, /i and g. The 
coefficients with L < n are governed by the renormalization group equation and can 
be determined recursively with the condition Ao,o = !> while the coefficients An,n(= 
An) need to be fixed by perturbative calculations and are called "non-logarithmic 
corrections." 



m 



n,L 



In 
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If we set the scale /i equal to m, which is accomplished at the specific point 

= m^{fi = M), (2-2) 
then all the logarithmic terms vanish, and we obtain 

m\M) (1 + A,g\M) + A2g\M) + ■■■), (2-3) 



where we denote the running mass and coupling constant at the scale fj, as m?'{^) and 
g'^ifi), respectively. Although both g'^(^A) and m^(7W) should be determined by the 
full renormalization group equation (RGE), we can approximate them by replacing 
them with gf_i^^p and rnf_^^^p, which are solutions to the RGE at one-loop order, 
namely 



m 



l-loop 



l + bogHf^)H- 



^l-loop ' 



70 ) 
2bo 



2 

9l-loop 



l + bog\fi)H- 



'l-loop ' 



(2-4) 



where bo and 70 are the coefficients of the renormalization group functions (3{g'^) 



-2bog^ - 0{g^) and j^{g^ 



log 



0{g ). Thus we have 



m 



l-loop 



,1 + ^l-loop^l + 5'l-loop^2 + 



(2-5) 



According to Appendix [HI this approximation corresponds to the leading-logarithm 
approximation with perturbative non-log corrections. At this stage, we assert that 
the approximation is valid if the running coupling constant gf.i^^p remains small at 
H = M.. Thus, we can apply this approximation to asymptotically free theories, if 
Ai is large enough. Thus we should choose the initial value of the running mass 
m'^{A4) in order for these conditions to be satisfied. 

2.1. One-loop improved RG analysis 

Let us consider a pole mass at one-loop order in ordinary perturbation theory. 



Mi 



m (fi) 



l-7lg2(fl)ln{^]+A,g\fl) + 0{g 



m 



(2-6) 



where 70 must be such that this expression is consistent with the expression (j2-5p . and 
^1, which comes from a finite part in the renormalization prescription, is the non- 
logarithmic correction. Note that we have used a "mass-independent renormalization 
scheme", like the MS or MS scheme. 

The leading logarithm contributions can be included as explained above, and 
then the pole mass becomes 



- TO 



1 + Wln(- 



loop ' 



+ "^l-loop ("^)5Loop I 



m)Ai. 



(2-7) 
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It is convenient to rewrite this as 



Mi 



m? (/i) 



1 + 



A' 

boF 



where 



1 



+ In 



m 



l-loop 



Note that F satisfies the recursive relation 

In ~'o 
F = -\ni[bo9\fi)F]^'o 



A2 



MS 



(2-8) 



(2-9) 



(2-10) 



with the basic scale = /x^e ^os^, where = Ane "'ji^. Let us define the 



_ 70 



dimensionless parameter x = {b^g^) 26n„^2 / /i2 
mass and the above defining equation for F become 



om Then the expression for the pole 



X 



TO 
P2b0 



1 + 



A' 



and 



F = -^hiF + hix. 
2bo 



(2-11) 



(2-12) 



Note that this function has a logarithmic cut starting from x = in the complex 
x-plane. We choose this cut so that F is analytic for a; > 0. 

Now we have a renormalization group expression for the pole mass at one-loop 
order, and the next task is to apply the variational method around a massless theory 
in the manner reviewed in the introduction. 

The improved renormalization group analysis for the massless theory is realized 
by the substitution 



m 



(1 - A)m2 
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(2-13) 



in the parts that admit power series expansions of integer powers in A. Note that we 
have improved the pole mass expression by use of the renormalization group, and 
actually its recursive form is a result of this refinement. As a result, there are some 
terms in Mp that do not admit a Taylor expansion. Because the coupling constant 
appears in a nonperturbative fashion, it is not difficult to see that it is sufficient for 
the improvement to replace x with (1 — X)x and to leave the rest unchanged. Then 
the mass Mp comes to depend on the formal coupling A and is expanded in power 
series in A as 



m2 (A) = J2 «-A^ 



(2-14) 



m=0 
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Fig. 1. Contour path. 
We define the n-th order mass after setting A = 1 as 

n 

K) = E 



m=0 



dz fl 1 
- + ••• + 



2'iri V z 



71+1 



1 + 



A' 



(2-15) 



Note tliat tlie analyticity of F leads to the conclusion that [F{(1 — z)x)] '^''o and 

[F{{1 — z)x)] ^''0 cannot be singular at the origin of the complex z plane when 
X > 0, and therefore the first term in the integrand must vanish. Let us define 
u = n(l — z). Then we have 



M, 



in) - ^^MS 



du / u\ -("+!) 1 
2^i\~n) ^ 



1 + 



A' 

6nF 



(2-16) 



Here, the contour path L is taken around the cut on the negative real axis in the 
u-plane, as shown in Fig. [TJ 

The approximation becomes better as higher-order terms are taken in the per- 
turbative expansion in the formal coupling A. Thus it is natural to take the limit 
n — > CO, and we then obtain 



Z,1lb JP2bQ 



1 + 



bnF 



(2-17) 



From the defining relation, F can be expanded around x = (hence u = 0) and we 
have 



[Fin)]- 



IXU) 



1 + (xu 



,1/B 



+ 



B -2 
2B 



(xu 



(2-18) 



where B = 7o/(26o)- We find an approximate expression for the pole mass with the 
help of the formula 

1 _ f dt 
Tiz)' Jl2^i 



(2-19) 
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This expression is 



M2(x) l + BAi Ai x-V-B / B'^-lAi\ xV-B 



B bo bor{i + i) V boJr{i-i) 

We would hke to find a plateau with respect to the parameter x0 because the 
dynamically induced mass can be accurately approximated by the value on such a 
plateau. 

We can further refine the above result by introducing a scale-changing parameter 
a, which is defined as a scaling parameter that rescales Ji to a/x. Although the 
rescaling of a renormalization point does not affect the renormalization invariant 
quantities, such as the one-loop mass rn^ Yf^^^, the perturbative mass does depends 
on the parameter a. However, the pole mass must be independent of a, and thus it 
is natural to tune a in such a way that a plateau clearly emerges. Once we introduce 
the parameter a, the defining equation for F becomes 



F = lnf^)-^lnF. (2-21) 
V a / 25o 



Finally, we arrive at 



M% f du „ ^-Jo. 

P - ^ I ^up 2bo 



Ay 2m 



(2-22) 



Thus, if we know bo, 70 and Ai, we can calculate Mp/A^^-^ as a function of x 
and a. We then search for a plateau of this function with respect to x by varying 
the parameter a and take the value on the plateau as the approximated value of 

§3. Improved RG analysis for the nonlinear sigma model 



Having reviewed the basic strategy, we move on to the 0{N) nonlinear sigma 
model in two dimensions as a preliminary example. The 0{N) invariant Lagrangian 
density with an external field is 

£ = ^{d^S-d^'S-2h-S) , (3-1) 
2g^ 

where S is the sigma model field of magnitude 1 (i.e. S ■ S = 1), and h is an 
A^-dimensional external vector field. It is convenient to choose h = (/i, 0, . . . , 0) by 
0{N) symmetry and decompose S as S = {a = — 7r2,7r*), where vr* are the 
components of an (N — l)-dimensional vector field with the condition vr^ < 1. Then 



Recall that x corresponds to the mass parameter to be tuned. 
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the action can be rewritten as 



1 

1 

2^ 



l-7r2 

d'x { {d^TT^r + (tt • d^nf + 



2ha 



2h + hn^ + -{TT^f + ■ ■ ■ } , (3-2) 



where the external magnetic field h plays the role of the infrared regularization to 
give a mass to the vr fields. We regard this as a massive counterpart of the original 
0{N) nonlinear sigma model without the external field. 

The IPI two-point function 77*-' of two vr fields at one-loop order is 



+ 



N 



-h 



/i-^/2(47r)^/2-ir(e/2), 



where d = 2 — e. It follows that the physical pole mass is 

^2 



= h 



47r 4 



(3-3) 



(3-4) 



It is easy to see that there is no non-logarithmic perturbative correction at one-loop 
order. 

We also obtain the renormalization functions 



and thus we have 



N 



2tt 



7n 



N 



An 



70 



+ 1 

Svr 



(3-5) 



(3-6) 



Substituting these into the pole mass formula (j2-22p . we find the plateau (in the case 
= 3) depicted in Fig. [21 in which we have plotted Mp/ with respect to h with 
various values of a. We recognise that a plateau emerges when we set o = 1, and the 
value on the plateau is Mp/Aj^ = 1. This result is indeed identical to the exact 
value in the large- A'^ limit, M^/yl^ = 1. Note that also in the Gross-Neveu model, 
the one-loop result is exact in the large- A'^ case.^^^ However, it is rather far from the 
exact value for = 3, M/Aj^g = 8/e ~ 2.940 It can be seen that the deviation 
from the basic scale A depends on the value of the non-logarithmic perturbative 
correction. It has been previously shown that the two- loop results for the Gross- 
Neveu model exhibit good agreement with the exact results even for smaller A^.^''^ 
Thus, we believe that if we proceed beyond one-loop order, we will obtain a more 
accurate result in this model as well. 



*' The expression for the exact mass gap evaluated in Ref. 29) is 



M = 



(3-7) 
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Fig. 2. The plateau in the nonhnear sigma model. The plateau emerges at Mp = 1 when 
a — 1. The horizontal and the vertical axes denote h and M^/A^, respectively. 



Let US comment on the spontaneous symmetry breaking here. Usually a non- 
zero expectation value of a mean field after the appropriate subtraction indicates 
that symmetry is broken spontaneously. However in the non-linear sigma model in 
two dimensions, there never occurs a spontaneous symmetry breaking because of its 
dimensionality. Thus a non-zero value of a mean field does not imply SSB. Rather, it 
merely indicates that the starting point of the perturbation theory does not respect 
the symmetry. In the non-linear sigma model, the generated mass term is indeed 
0{N) invariant, and thus SSB does not occur. If the generated mass term breaks a 
symmetry, then one can conclude that SSB takes place and use the generated mass 
term as an order parameter, as in the Gross-Neveu model, where a fermion mass 
term breaks the chiral symmetry. 

§4. Improved RG analysis for Yang-Mills theory 

Let us start with a comment on what we are going to calculate using the improved 
RG analysis in pure Yang-Mills theory. The improved RG analysis presented here 
is basically a method to calculate a mass term in a theory. In Yang-Mills theory, 
the mass term for the gluon is known to be strictly forbidden by gauge symmetry. 
Indeed, it is shown in Ref. 25) that there is no mass dimension two BRST invariant 
local operator which can be a source of mass generation. In a specific gauge, however, 
there can exist a local mass dimension two operator which may be a source of mass 
generation, though this "mass" does not have any gauge-invariant meaning and is 
not a physical mass. Below we consider one of these operators, TiA^A^, in the 
Landau gauge. This is gauge dependent and becomes nonlocal in a general gauge. 
However, as long as we start with the gauge-fixed action in the Landau gauge, in 
which there is no classical gauge invariance, we can think about the condensation 
of this operator, and through the four-point interaction of gluons, a mass term for 
the gluon may be generated. In the gauge-fixed action, BRST invariance is an 
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important symmetry, and therefore it has to be maintained throughout calculation. 
This operator, with a space-time integral, is indeed BRST invariant in the Landau 
gauge. In the following section, wc consider the Curci-Ferrari model as our massive 
gauge theory counterpart. Although it seems at first sight that the mass term will 
undermine the unitarity of the original Yang-Mills theory, our mass term is fictitious. 
In the prescription of the improved RG method, we introduce a counterterm which 
cancels the effect of the mass term. Therefore, we expect that pathologies of the 
massive gauge theory will eventually disappear. 

As we will see, the mass term has a nonperturbative nature; that is, it is propor- 
tional to the scale given by the dimensional transmutation. A proportional coefficient 
would be of order 1. This term might be corrected by perturbative quantum effects, 
like the instanton action. However, due to the lack of a clear understanding regard- 
ing the mass generation and condensation of this operator, it is difficult to calculate 
corrections to determine the exact amount of the generated mass. For this reason, 
wc restrict our investigation to the leading order, and the formulation of a precise 
argument for these coefficients is postponed to a future study. 

The Lagrangian of the Curci-Ferrari model that is our counterpart of U{N) 
Yang-Mills theory is 

C = -^{d,Ai-d.A-f 

- ^{d^A-n^ - ^gr'^'^A^.AUd^A^c - d^A^c) - f"^'' A^^A^.A^^^A^ 



,2 



8 

+ ^n?AlA''^' - ml^^c" , (4-1) 

where ^ is the gauge parameter and /? is an extra parameter characteristic of the 
model. Here, the ghost mass is related to the gluon mass as rn^^ = ^m?. The 
condition that the model possesses the BRST symmetry is given by 

^Pni^ = 0. (4-2) 



We carry out the calculation by setting /3 = 1. The BRST invariance will be manifest 
in the last stage of the calculation, where we could regard the original mass as going 
to zero as a result of the subtraction or when we go to the Landau gauge, ^ = 0. 

To find the perturbative expression for the pole mass of the theory, it is sufficient 
to calculate the one-particle irreducible (IPI) two-point function. We write the IPI 
two-point function as follows: 

in^^ = idab [n.m^v'"' + TTsg'^'^l , (4-3) 
where Af^" = rj'^" — jf^ ^ and tti and 7r2 should be calculated perturbatively. 
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Then, the full propagator iE^'^ is written in the form 



-i/jl-TTi) 

^ l-7r2 



l-7r2 



+ 



-i/(l-7ri) 



Thus, the pole mass can be read off as 



Mj, 



^2 1+^ = m\l + TTl + 7,2) + 0(/). 



(4-4) 



(4-5) 



I - TT2 

The IPX two-point function is calculated up to one-loop order in a standard way 



as 



in: 



ab W 



where e ^ = e ^ 



dx / dy 



dz 



AO) 



(47r)2 

- In Aim^T^^^NmO + q^A'^'^NAo) + 2im^rj^''Nmi + q^A^^NAi)] , (4-6) 
7 + ln(47r), 

Z\ = [1 - (1 - e) {zx + (1 - z)y}] - z(l - (4-7) 

and the explicit expressions for the coefficient functions N^oix, y, z; ^), N^iix, y, z, m?, q^'-,£,)i 
NAo{x,y, z;(,) and NAi{x,y, z,m? , q'^; S^) are listed in Appendix [Al 

The constant Ai is given by the finite part of tti -|- 7r2, which is found from (j4-6|) 
to be 



lN_ 

(4^)2 



313 
"36" 



11^: 



•vr 



7 (4^ _ 1)3/2 

— t H arctan 

12^ 12 



1 



+ 



Qi- 1 

24 



In the Landau gauge, we have ^ ^ 0, and Ai becomes 



Al 



(4^) 



313 ll^/37r 
"36" 8 



The renormalization functions are 

^2 



0(9' 



7r) 



35 - 3e^_ 5^ 



12 



(4vr) 



+ 0{g' 



and thus 



11 iV 
"6"(4^ 



70 



35 - 3C iV 
12 (47r)2 



In^ . 
(4-8) 

(4-9) 

(4-10) 
(4-11) 



Using the above results, we apply the improved renormalization group method 
explained above to Yang-Mills theory and are able to calculate the approximated 
pole mass of the gluon. Note that at one-loop order, the result is independent of 
N. We plot Mp/yl^^ with respect to x, which corresponds to the provisional mass 
771,2, £qj, various values of a in Figure [31 Figure [3] shows that a plateau emerges near 
12_ 

MS 



M2/vl2 ^0.66. 
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X 



Fig. 3. The plateau in the Landau gauge. The vertical axis represents Mp/A^j-g, evaluated using 
(I2-22|) with the data (|4-9p and (|4-lip . The horizontal axis represents x proportional to m^. We 
have plotted the graphs for various values of a from a = 0.6565 to a = 0.6635. The plateau 
appears near Mp/A^^ ~ 0.66. 



§5. Conclusion 

In this paper, we have presented one method for evaluating the operator A = 
I (volume) ~^(min(/ / d^x{A^)'^). In the Landau gauge, the calculation boils down to 
evaluating {A^A'^) , which may be connected to the mass pole of the gluon propagator 
in this gauge. We have thus applied the improved renormalization group method to 
pure Yang-Mills theory at one-loop order and obtained the result that a nonpertur- 
batively generated pole mass of the gluon, Mp, emerges as 

M2/yl2_ ~ 0.66, where 

-^Ws scale. Some work has been done to calculate the mass of the gluon, 

for example, in a lattice calculation. ^^-^ Our result is much smaller than those found 
in the previous work. However, we cannot compare the previous result with our 
result directly, because that result is not the pole mass itself but, rather, the vacuum 
expectation value of the A^A^ operator calculated through the operator product ex- 
pansion. It would be interesting to perform further calculations, two loops or higher, 
and observe if higher-order calculations might give larger gluon masses. Calculations 
of the renormalization functions of the Curci-Ferrari model have been carried out 
to three-loop order. ^''^ Those calculations enable us to improve the renormalization 
group analysis, but we think that perturbative corrections play an important role in 
producing a nonperturbative mass. Again, we would like to emphasize that though 
we have done the calculation in the Landau gauge, the result should have a physical 
meaning, because the original Z\ is a gauge-invariant quantity. 

Here, let us comment on a possible relationship between mass and confinement. 
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If the pole mass of the gluon were infinitely large, we could interpret it as a signal 
of confinement, since such a large mass would indicate that the gluon does not 
propagate as a physical mode. This is an interesting viewpoint^ and we expected 
to find a large gluon mass. Although one-loop order we have not obtained such a 
large mass, the gluon mass might become larger and larger as the order of calculation 
increases. 

Yet another perspective concerning confinement may exist. In the operator for- 
malism of Yang-Mills theory, there is a well-known criterion for confinement, the 
so-called Kugo-Ojima condition, ^'^^ which is a sufficient condition for color confine- 
ment. The condition mentioned here concerns a BRST transformation property of 
anti-ghost fields, and the criterion asserts that when this condition is satisfied, one 
can prove that no color singlet state belongs to the BRST cohomology, and hence it 
cannot be observed as a physical state. This condition has been interpreted as involv- 
ing the infrared properties of gluon and ghost propagators in the Landau gauge, 
and these properties have been tested with lattice simulations,^^-* analytically^'^^ and 
by employing the Schwinger-Dyson equations. Because the infrared behavior of 
the propagator is involved with the existence of the nonperturbatively generated 



mass, it would be interesting to study it with our method!**^! This viewpoint also 
provides an interesting interpretation of the remark above, that is, that the confine- 
ment via the Kugo-Ojima mechanism requires a finite mass for gluon propagators, 
and thus it may be consistent with our result. 

Another issue concerns the ghost mass generation. In the literature, generation 
of the ghost mass has been studied intensively, and it might also be possible to 
evaluate it using the prescription presented in this paper. It would be interesting to 
calculate the pole mass of the ghost propagator. 

We have also applied this method to the 0{N) non-linear sigma model in two 
dimensions. The exact mass gap of this model is known, and we can compare 
it with our result in order to check the effectiveness of this method. However, at 
one-loop order, the calculation is consistent with the exact result only in the large- A'^ 
limit, where we have Mp = ^jg^. A similar fact was observed in the Gross-Neveu 
model^^^ at one-loop order. In Ref. 17), a two-loop calculation is performed and 
a result in good agreement with the exact result, even for small values of A, is 
obtained. Considering this fact, it would be interesting to calculate the next order 
in the non-linear sigma model to make sure of the validity of this method. 
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Appendix A 

Coefficient Functions 



Below we give explicit expressions for the coefficient functions A^^o , ^Ao , ^mi 
and NaV- 



N^o =-z{l - z) (16(1 - (,)zx + 16(1 - e)(l -z)y-b + 9^ , (A-1) 

Nm) =z{l - z)(80z(l -z)- 3(1 + 0) , (A-2) 
A^i = [{-19(1 - ifz{l - z)((l - z)y + zxf 

-^(63^ - 181)(1 - 0^z{l - z){{l - z)y + zxf 
-^(1 - 0(2^' - 35^ + 45)^(1 - z){{l - z)y + zx) 

-^(l-0(2e-5Ml-^)}(m')' 

+ {(1 - 0^(-46z(l -z) + 3)z(l - z){(l - z)y + zxf 

X (-174^2(1 - z)+ 446z(l - z) - 3C - 21)z(l - z)((l - z)y + zx) 
+\z{l - z)(-24^2^(l -z)- 76z(l -z) + 156^z(l - z) - 9^ - 3)| rn^q^ 
+ {(1 - 6(-26z(l -z) + Z)z\l - zf{{l - z)y + zx) 
+-z^{l - z)^(-72^z(l -z)+ 32z(l -z)-l + Q^) \ {q 



.2\2 



IA\ (A.3) 



Nay 



-(-116z(l - z) + 7C + 10)z(l - z){{\ - z)y + zxf 

+(1 - 0(-22ez(l - z) + 122z(l -z)^i^ -li- 11) 
X2;(l - z)((l - z)y + zx) 

+^z(l - z){Vl^z{\ -z)- 258z(l - z) 
+ 86^^(1 -z)- 3^2 + 15^ + 24) I {m^f 

^(1 - ifz{\ - z){\z{\ - z){272z{l -z)- 33) 
-28^^(1 - z) - ^2 _ 4^ _ 3)(^i _ ^ 
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+-z{l - z){4z{l - z)(281z(l - z)- 34) + 4^z(l - z)(33z(l - z) + 10) 



+ 8^2^(1 - z) + 3^2 + ee + 3)| m 

+ ^z\l-zf{8^z{l-z)+e + A^ 

+3 - 8z(l - z)(76z(l - z) - ll))(g^)^] /A^. 

Appendix B 

Renormalization Group Analysis for the Pole Mass 



(A-4) 



In this appendix, we explain the fact that a pole mass can be written in terms of 
a logarithm of m? / fj?, following Ref. 28). We also comment on a leading-logarithm 
approximation. 

Here we concentrate on a theory that has two dimensional parameters, and 
//. The physical mass that appears at a pole of the propagator is 



m 



n=0 



In 



(B-1) 



where the quantities fn are undetermined functions of w? j [i? and are calculated by 
perturbation theory, and g is the coupling constant of the theory. Note that here 
we have applied a mass-independent renormalization (MIR) scheme, like minimal 
subtraction MS or MS. Only when an MIR scheme is applied does the pole mass 
have a form like ()B-ip . 

By definition, the pole mass Mp is independent of [i\ that is, it is renormalization 
group invariant. Then, Mp satisfies the renormalization group equation (RGE) 



+ /3(5')^ - 7m2(5)"T'^^-^ 



mA = 0, 



where 



Substituting (lB4]l into this RGE, we obtain 



/i dm? 
m^ dfi 



(B-2) 



(B-3) 



In 



Eh 



.n'<n 



d 



,2m 



Integrating with respect to In/i, we obtain 
d 



fn'9 



2n' 



cooff. of g 



fn 



E 

.n'<n 

-|- const. , 



7m2m 



(B.4) 



coefT. of g 



(B-5) 
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where const, is independent of fi and, thus, m^. We have the "boundary condition" 
/o = 1 of this integral equation, and hence 



/i = const. 7o In 



m 



where we have defined 



Pig) = -26o/ - 26i/ + 



705^ + 7i/ + 



(B-6) 



(B-7) 



In this manner, we can determine all the functions /„ recursively by using of infor- 
mation provided by the RGE, up to constant terms. Thus we obtain 



Mi 



m 



71=0 L=0 



n—L 



9 



2n 



(B-8) 



where the coefficients An^L with L < n are determined through the above procedure, 
while An^ni= An) are unknown constants at this stage. These constants are to 
be calculated within the perturbation theory, and we call them "non-logarithmic 
corrections" . 

An interesting fact is that the L = contribution in the series (jB-Sp . called 
a leading logarithm approximation, satisfies the renormalization group equation at 
one-loop order, that is. 



d_ 

djjL 



, 2 ^ 2 2^ 

2"o5 -K-n + log m 



0. 



Therefore we can solve it in a standard way, and we obtain 

.2^ 



LLn2 



(M^L) 



m 



oo 
n=0 



A^ 



n.O 



g 



2n 



(B-9) 



(B-lO) 



This indicates that once we know and 70 at one-loop order, we can calculate the 
leading-logarithm contribution to the pole mass at all orders in g"^ . 
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